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Abstract 

We study the equations for the quasi-momenta which characterize the wave-functions in 
the Bethe ansatz for the XXX-Heisenberg model. We show in a simple analytical fashion, 
that the usual "string hypothesis" incorrectly predicts the number of real solutions and the 



number of complex solutions for N > 21 in the sector with two spins flipped, confirming 
the work of Essler et al. Two complex pair solutions drop out and form two additional 
real pair solutions. We also introduce a new set of variables which allows the equations to 
be written as a single polynomial equation in one variable. We consider in some detail the 
case of three spins flipped. 

The XXX-Heisenberg model appears to be fundamental to many of the recent interest- 
ing developments in modern quantum field theory and mathematical physics. Integrable 
models, conformal field theories and quantum groups are seen to arise from various limits, 
extensions and generalizations 1 . It is, in fact, the original model considered by Bethe 2 , 
which gave rise to the celebrated "Bethe ansatz" solutions. Recently, there appeared an 
article 3 , studying the Bethe ansatz equations in the two particle sector of the spin one- half 
model. It was shown that, as the number of lattice sites, N, increases past 21, two new real 
pair solutions appear. These solutions do not fit into the conventional scheme of classifi- 
cation of the solutions to the Bethe ansatz equations which causes some anxiety as to the 
verity of the completeness of the full set of 577(2) extended Bethe ansatz states. The full 
set of states, just by counting, should number 2^. This, however, is not a problem. As the 
real pair solutions appear, complex pair solutions disappear, conserving the total number 
of states. We analyze the Bethe ansatz equations analytically and confirm their result; two 
complex pair solutions disappear simultaneously giving two new real pair solutions. 

The spin \ XXX-Heisenberg model with N spins corresponds to the Hamiltonian 




(1) 



with <tat_|_i = o \ , and the <7j are just the Pauli matrices for each i. The Bethe ansatz in 
the M particle sector corresponds to eigenfunctions of the form 

l*)= ^i r --,XM)<T~<x~<x~ ••■<7x M \°) ( 2 ) 

x 1 <x 2 <---<x M 

with 

i)(x 1 ,---,x M )= Yl e PU)>p(l) ' ( 3 ) 

PeS M 

|0) is the ferromagnetic state with all spins up, P is a permutation on M objects and dis- 
satisfies 

2cot(^)=cot(|)-cot(|). (4) 

k\ ■ ■ • ku £ [0, 27r] are the quasi-momenta or spectral parameters. The wave-function is 
symmetric under interchange of any two coordinates Xj. 

Periodic boundary conditions 

ip(xi, ■ ■ -,x M -i, N + 1) = ip(l,xi, ■ ■ ■ ,x M -i) (5) 
imply the M coupled equations 

M 

= JJ e ^,« j = l,2,-..,M. (6) 



e ikjN 



1*3 



The usual set of variables considered are Aj = cot(-^) which give the equivalent form for 



the equations 

(e(A,)f = jW^M J = 1,2,---,M, (7) 



i=i 

1*3 



with 



, . . A + i . . 

e(A) = —.. (8) 
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The energy of the state is then given by 



M -2T 

E(A 1 ,-.-,A M ) = J2 J T T J- (9) 

j=l 3 



In the large N limit, the "string hypothesis" 4 states that, for fixed M, any solution 
Ai, ■ • • , Am consists of "strings" of the form 

= AS + i(n + 1 - 2j) + o(e- SN ) j = 1, • • • , n (10) 

where n > 1 gives the length of the string, a labels strings of a given length, j specifies the 
imaginary part of A and 5 > 0. With such a hypothesis, the Bethe ansatz equations only 
involve the real parts, A™. The solutions are parametrized by (half-odd) integer numbers 
I£ for N — M n (even) odd, where M n is the total number of strings of length n. Clearly 
S^Li n M n = M. It is generally believed that there is a 1 — 1 correspondence between 
solutions of the Bethe ansatz equations and sets of independent, non-repeating integers 
^, (7™ 7^ Ip for a ^ /?, i.e., no two strings of the same length contain the same integer, 
within one solution set of Ai • • ■ Am) 4 . 

The Bethe ansatz equations for M = 2 are 



(e(Ai)) =e( ) 

( i a / A 2 - Ai 
(e(A 2 )) =e( ). 



(11) 



This set of equations has the following symmetries, Aj — > A*, Aj — > — Aj and Ai <-> A 2 . 
Now 

Al _ A Ai^+i _ Al + Z _ ( A 2 _,) 

1 2 } Ai-i-CAa + t)' 1 J 
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Replacing for Aj with Aj = l -^x^T gi yes 

rr Ai-A 2 ,_ / e(A 1 )e(A 2 )-2e(A 1 ) + l \ 
1 2 ] Ve(Ai)e(A 2 )-2e(A 2 ) + i;- 1 ^ J 

We choose to work with the variables Xj = e(Aj). In terms of these variables, real solutions 
(in the Aj) are mapped to the unit circle while complex conjugate pairs are mapped to 
complex pairs which are related by, (z, j^). Then we get 

-N ( XiX 2 — 2Xx + l\ N ( X\X 2 — 2X 2 + 1 



*■ =- u;*;-2*;+i j and ^- Ua-^+j - (i4a '" ) 

Multiplying these equations together gives, 

(X^)^ = 1, hence XiX 2 = w (15) 
where uj is an Xth root of unity. Thus replacing X 2 = in (14a), yields 



Assuming Xi 7^ ^j-, or 0, we get 



Xf^ 1} ((w + l)Xi - 2co) + to - 2Xi + 1 = (17) 



le. 



(oj + l)Xf - 2uX[ N 1} - 2Xi + (u + 1) = 0. (18) 

We will be interested in the real roots in terms of Ai. These are mapped to roots on the 
unit circle in terms of X\. Furthermore 



2u 



(00+ + 1) 2 + 2Rea; 1 + Reu; 
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> 1 (19) 



for (jj 7^ 1. Thus the denominator by which we multiplied may only vanish for u = 1 and 
X\ = 1, which gives X 2 = 1. The corresponding wave function vanishes identically as the 
Bethe ansatz wave functions respect the Pauli principle. We want the roots of equation 
(18), for each choice of u. This gives N polynomials, each of order N, which is a total of 
iV 2 roots. X\ = y/u, however, is always a solution with X 2 = y/u. These wave functions 
also vanish because of the Pauli principle. Furthermore, if X is a root, then ^ is also a 
root. Each one gives the same Bethe ansatz wave function, only the roles of X\ and X<i 
are exchanged. Therefore, we get a total of n2 ~ n = (^) different wave functions. This is 
exactly the dimension of the subspace of states with two spins flipped. There are exactly 
( ^ ) independent ways of flipping two spins among N. 

The equation is cast in a more symmetric form with the replacement X\ — > y/ujX\ and 
correspondingly Y\ — > y/uY\. These satisfy the symmetric relation Y\ = We first 



consider the case N odd, where it is always possible to take y/uj N = — 1. This means with 
oo = e % ^~ , for m odd we take y/u = e l ~ but for m even we must take y/u = e~ % n . 
The set of {y/u>} f° r m °dd are exactly the inverses (complex conjugates) of the set for m 
even, except for u = 1, which is not so paired. Then taking into account y/u N = — 1, we 
get 

= (u + 1)V^ N X™ - 2uy/^ N ~ 1 X i 1 N ~ 1) - 2y/ZjX x + (u + 1) 
= -{u + 1)X? + 2y/^x[ N ~ 1} - 2y[ZjX X + (u + 1) (20) 

= -yfr ({yfa + -^=)X? - 2Xf" 1 + 2X X - + -L) 



which yields 



cos( d -)X» - X[ N ~V + Xl - cosi -) = 0, (21) 



with the definition y/cJ = e 1 ^). The symmetries of equation (11) now translate into 
(X 1 -> -> w), (X x -> -i- jU ; -> w*) and (X t -> ^-,cj -> w). 



The roots of Equation (21) come in pairs, if 2 is a root, so is |. Furthermore, equation 
(21) has real coefficients, thus complex solutions come in complex conjugate pairs. The 
natural grouping of complex solutions is actually in quartets, (z, z*,^, pr). Real solutions 
in terms of the A variables are mapped to solutions on the unit circle, thus the quartets 
degenerate into pairs for these. Complex solutions in terms of the A variables are then 
mapped off the unit circle, and should generally come in quartets. Therefore, if there 
is truly only one complex pair in terms of the A variables, it must be mapped to a real 
pair, r, £ in terms of the X variables. X± = r fixes the solution for the corresponding 
^ = Cv^+i = ~2 pqrf + ^T^TT- ^ would be interesting to find true quartet solutions. 
Below we only find pairs. Equation (21) is the same for ^Ju and for Although the 
equations, and consequently the solutions, are the same in the tilde variables, the solutions 
in terms of the original X variables are of course different, giving rise to different Bethe 
ansatz wave functions. 

Now we look for solutions on the unit circle, X\ = e ia . Equation (21) becomes 



cos Ae l7Va - e l(Ar - 1)a + e ia - cos A = 0, (22) 



which simplifies magically to 



iNa Q iNa iNa ■ i(N -2)a i(N -2)a -j( N -2)a 

= e 2 cos(-)(e 2 -e 2 ) — e e 2 (e 2 -e 2 ) 



iNa 



e 2 2z(cos(-)sin( — )-sinO 



(23) 



7 



Thus 



A sin(^) 
cos (o) = ■ 7lk\ 



(24) 



This equation is easily studied graphically, see figure (1). a G [0,27r]. The right hand 
side has zeros at a = for k = 1, 2, 3, • • • , N — 3 and poles at a = for k = 

1, 2, 3, • • • , N — 1. The value at a = is but at a = rc it is —1. We get a root for each 

. , (N-2)a s 

intersection of a horizontal line with y = cos(|) with the curve y = S1 "- n( - Na^ • 



N=ll 




Figure 1: /(a) 



sin( 



QV-2)c 



sin(^) ' 



JV = 11 



By counting the intersections it is evident from figure (1), that there is a root for each 
zero. There are always iV — 3 such roots. An additional pair of real roots can appear if 

(25) 

i\ 
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cos(-) > — — — . 
V - N 



giving N — 1 real roots, but X\ = 1, [X\ = ^/uj), is also a root of the original equation, 
giving totally A" roots. This exhausts all roots of the polynomial and there are no complex 
conjugate pairs (for A). Thus we find 

/ ^critical \ AT 2 /<->/•• \ 

cos( — - — ) = — ^— . (26) 

For large N, this will have solutions for | near or 2ir, but because of reflection symmetry 
about a = 7r we need only search near a = 0. Assuming 9 is small and making an expansion 
in 6 and we get 

1 1 /^critical n2 , _ -i _ ^ , 



giving 



4 

^critical = 7= • (28) 

viv 



The values taken by | (the condition that ^/uJ^ = — 1 must be satisfied) give 



TtTYl 

- = — for m=l,3,5,-..,iV-2, (29) 



or 



tt 777, 

- = 7r-— for m = 2,4,6,...,JV-l, (30) 



AT 

(A" is odd). Hence, for the odd series in to, we get 



2 /— 

m<-7iV (31) 



7T 



and for the even series in to, 



2 y— 

|A^-to|<-VA7. (32) 



7T 



The number of new solutions behaves like v A" 



The first N for which we have a new solution, (to equation (25) actually), is N = 23, when 
m = 3 is allowed, next at iV = 63, m = 5 is allowed and so on. (The solution at m = 1 is 
accounted for in the string hypothesis set of solutions. It corresponds to the action of the 
lowering operator applied to the corresponding state in the M = 1 sector, in the context 
of the 577(2) extended Bethe ansatz.) 

Analyzing in more detail the complex pair of roots, we insert X\ = r into equation (21), 
giving 

r ( r N-2 _ j\ 
co <)= >_! ' =/(*•)• (33) 



N=ll 
1- 




Figure 2: f(r) = r(r r "" 1 " 1) , N = 11 

It is easy to see that f(r) is monotone decreasing for r > 1. Therefore, it is monotone 
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increasing for < r < 1, as /(-) = /(r). The maximum occurs at r = 1, /(l) = 



N-2 
N • 



Evidently there are two roots for each < cos(|) < ^r^, as we can see from figure (2). 
This is exactly the same condition as equation (25), as expected. We see that for | near 
Tj- the solution (r, ^) is (oo,0), which implies A « —2 r r f^ ± i, as predicted by the string 
hypothesis. For | near zero, however, the imaginary part becomes arbitrarily small and 
the string hypothesis is grossly violated, finally to the extreme that the complex solutions 
drop out all together. 



For r < 0, 



\ r \N-2 , 1 

/(r) = - |r|( w ] (34) 



It is again easy to see that f(r) decreases monotonically from to —1 as r varies from — oo 
to —1 after which it rises monotonically to 0, since always /(£) = f(r). Thus, we always 
get exactly one complex pair for | e (f, 4r), as we can again see from inspection of figure 
(2). 

Finally we note that actually the complex pair solutions come in pairs. This is because of 
the pairing of ^Ju for a particular m odd with -4=^ for the corresponding m even. These 
yield identical equations (21) and (33). The corresonding pairs of complex pair solutions 
are, however, not the quartets to which we had referred earlier. Interestingly enough 
though, clearly both pairs drop out simultaneously when the critical condition equation 
(25) is staisfied. Thus actually two complex pair solutions of the "string hypothesis" drop 
out simultaneously and become two real pair solutions. 

We present below in somewhat less detail, the case when iV even. First of all, we have 
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(y/uj) N = (— l) m , for either choice for the square root, y/u = ie 1 ™ . Thus, for the case m 
odd, we get the condition as above, 

A sin(^^) 

cos U = ■ J a , - > ( 35 ) 
2 sir^^Tp) 



but here N = 2n. Thus, 



6 sin((ra - l)a) 

cos(-) = r^^. (36 

v 2 ; sin(na) V ; 



The R.H.S. is a function which has zeros at a = (Jf-i) ' ^ = 1? 2, 3, • • • , 2n — 3, excluding 
k = n — 1, i.e. 2n — 4 zeros, and poles at a = ■^ L , A; = 1, 2, • • • , 2n — 1, excluding k = n. 
The value at a = is The situation is as before except at a = 7r, the value is — IL ^, 
not —1, and the function "turns over". The function is reflection symmetric about a = it. 
We therefore obtain for each u> at least 2n — 4 real roots (the number of zeros) , but we get 
two additional real roots if 

|cos(|)|>^, (37) 

making a total of 2n — 2 = iV — 2. This condition can be satisfied for | near zero or near 
7r. We also of course get the reflections of these about ir hence we can restrict | G (0, 7r). 
AtA^ = 22,m = 3is allowed, and then at iV = 62, m = 5 is allowed. The condition 
becomes m < -\/N and liV — ml < -\/~N for large A 7 ". Thus for the case even, but 
m odd there can be two solutions which violate the string hypothesis correpondence. We 
also always have solutions of the original polynomial corresponding to X\ = ±1, (N even) 
which gives totally A^ roots, in which case, there are no complex roots for values of m 
satisfying equation (37). We notice that equation (37) is satisfied simultaneously by two 
values, TOi,m2 such that mi + m2 = N. (N.B. mi are odd, Af is even.) Hence we again 
lose two complex pair solutions, giving two additional real pair solutions. 

12 



For m even, the equation for Xi becomes 

= (u + l)yfa N X? ~ 2uV^ N ~ 1 X { 1 N ~ 1) - 2^X X + (u + 1) 
= (u + - 2 V ^X { 1 N ~ 1) - 2v^X x + {lu + 1) (38) 

= v4J f (v^+ -^)^i V - 2X 1 JV " 1 - 2X X + (v^+ -L) 

Replacing Xi = e ia yields 



cos Ae l7Va - e*^" 1 ^ - e ia + cos A = 0, (39) 



(with y/uj = e % 2 ) which simplifies to 

tNa ,9 ,, iNa -iNa ■ i(N-2)a i(N-2)a -i(N-2)a 

= e 2 cos(-)(e 2 +e 2 ) — e e 2 (e 2 +e 2 ) 

iN ex. 



,6. x iVa N ,(N-2)a. 
e 2 2 ( cos(-)cos(— ) -cos( ) 



(40) 



yielding 



.6 = cos(_^^ 
2 ] cos(^) 



= ~TWorT^ • (41) 



This equation is also easily studied graphically, see figure (3). The R.H.S. has zeros at 
a = for k = 1? 3 5? . . . ^ 2Ar _ 5 and poleg a t a = *2i for jfe = 1, 3, 5, • • • , 2JV - 1. We 

can check that at a = the function is 1 and at a = 7r it is — 1. Graphically we get a root 

/ (N-2)a \ 

for each intersection of a horizontal line y = cos(|) with the curve y = c , . 

^ COS ^ 2 ) 

Counting the intersections it is evident from figure (3) , that there is a root for each zero. 
There are TV — 2 zeros, thus there are always only N — 2 real roots, X\ = ±1 are not 
roots in this case, in fact it is easy to see that cos and cos ^ ( JV ~ 2 ) Q j never vanish 

simultaneously. Hence there are always two complex roots and no new real roots which 
violate the string hypothesis correspondence. 
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We summarize in the following way. For N sufficiently large, we have 2 y^y/Nj additional 
real pair solutions with respect to the prediction of the string hypothesis. For N odd, they 
are distributed evenly between m odd and m even, however, they occur only near X\ = 1. 
For iV even, they occur only for m odd, however, now they are distributed evenly near 
Xi = ±1. We give simple, exact expressions for the critical values N, equations (25) and 
(37), which do not appear in Reference (3). 

Consider now the sector with M spins flipped. Here the equations (7), rewritten in terms 
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of the variables Xi = e(Aj) are 



Multiplying the M equations together gives 

(X 1 X 2 ---X M ) N = 1. (43) 

Thus 

X M = - (44) 

AxA 2 • • • Am-1 

with cu^ = 1. Multiplying through with the denominator we get the coupled polynomial 
system 

M M 

X l \[(X l X l -2X l + l) + (-l) M \[(X l X l -2X l + l) = i = l,...,M. (45) 

(=i i=i 

Using equation (44) we can eliminate Xm from the first M — 1 equations and multiplying 
through with the denominator yields the system of M — 1 polynomial equations 

M-l 



X t J[ (XiXt - 2X l + l)((Xi - 2)oj + X X X 2 ■ ■ ■ X M -i) 



i=i 

M-i ( 46 ) 

+ ("1) M II ( X * Xl - 2Xl + l ^ X ^ - ( 2X * - l ) X ^ • --Xm-i) = 

i=i 

i = 1,...,M- 1. 

Each equation is just a permutation of the variables in any other one. Removing an overall 
factor of Xi, the ith equation is of order N + M — 2 in Xj, but only quadratic in all the 
other variables. 
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(47) 



Specializing to M = 3 we find the system 

= X?{X X X* - 2X 2 + l)((Xi - 2)u + X X X 2 ) 

- (X X X 2 - 2X X + l)(X 1( v - 2Xl + X1X2) 
= X 2 N (X 1 X 2 - 2X 1 + 1)((X 2 - 2)u + X 1 X 2 ) 

- (X 1 X 2 - 2X 2 + l)(X 2 u - 2X\ + XiXz). 

Simplifying we get 

= X2(X? +1 - 2X? + 2X\ ~ Xi) + X 2 (u(X? +1 - 4Xf + AX» ~ x - X x ) 
+ (X? - 4Xf + 4Xi - 1)) + u{X? - 2Xf" 1 + 2Xi - 1) (48a, b) 

X x <-> X 2 . 

There are two ways to proceed, both give the same result up to trivial factors. The simplest 
and straightforward method is to solve the quadratic equation (48a) for X 2 

A(X 1 )Xl + B(X t )X 2 + C(X 1 ) = (49) 

i.e. 



Y _ -B(Xi) ± y/B 2 (X 1 ) - 4A(X 1 )C(XT) 

X2 ~ 2MX[) (50) 

as an algebraic function of X\. Then replace this in the corresponding equation (48b) for 

X! 

A(X 2 )Xf + B(X 2 )X 1 + C{X 2 ) = 0. (51) 

Expanding the powers of X 2 we can isolate the odd powers of the square root to one side 
of the equation, the other side then only involves rational polynomial expressions in X±. 
Finally squaring both sides and multiplying through by the denominator gives a single 
polynomial equation for X\. Generalizing from calculations with Mathematica for low N, 
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we obtain the form of the polynomial as (X x - l) Ar+4 X 1 JV - 2 A(X 1 ) JV + 1 P JV (X 1 ), where P N 
is of degree iV(iV — 1) + 1. It is not illuminating to present the explicit form of P/v for a 
few explicit value of iV, we do not have a closed form for it. The roots of the polynomial 
along with the equivalent of equation (44) should give rise to all solutions of the Bethe 
ansatz equations in this sector. 

The second way to proceed affords generalization to higher values of M, although this is 
not straightforward. We consider the quadratic equation in X 2 as giving the first step of 
a recurrence relation 

X 2 - B{Xl) X C{Xl) (52) 

x >--Mx7) X2 -Mx7y (52) 

Then if X 2 = a n X 2 + (3 n we have 

a n X 2 + (3 n = X 2 = X 2 X^~ X = X 2 (a n -iX 2 + fi n -i) 

= (a n -xX\ + fi n _ x X 2 ) 



B{X 1 ) v C(X 1 ) 
- («-i("^y x 2 - -j^) + (3n-lX 2 ) 



(53) 



«n^ _ 1 ( -B{X 1 ) A(X 1 )\fa n _ 1 \ _ M /a n _i , ( , , 



Thus 

(p n j~ A(Xi)\ -c(xo o ; yp n -ij ~ J " 1 ypn-! 

which is trivial to solve as 

A(X 1 ) 
CjXi) 
A(X 1 ) 

Then replacing in equation (51) for X 2 we obtain a linear equation in X 2 that is a rational 
polynomial expression in Xi and we solve for X 2 as a rational polynomial in X\. Replacing 
this back in equation (49) and multiplying through by the denominator again gives a 
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Z) =Mn ~ 2 -it) ) (^) 



single polynomial equation in X\. We obtain from Mathematica the same polynomial as 
the above, modulo some trivial factors, (Xi — 1) N+1 X( J A(Xi) N Pn(Xi). However, since 
we have no closed form expression for P/v, whose degree is of order iV 2 , it is actually 
quite unfeasible to go much beyond N = 20. We are presently engaged in numerical and 
analytical analyses of the roots of this polynomial to see if there are any new violations of 
the string hypothesis. 

We find it is a dramatic simplification to deal with even, coupled polynomial equations 
(46) than the original transcendental equations. We hope to extend our analysis to other 
models where the Bethe ansatz has proven useful. 
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